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Abstract 

In this paper, the analog of Maxwell electromagnetism for hydrodynamic 
turbulence, the metafluid dynamics, is extended in order to reformulate the 
metafluid dynamics as a gauge field theory. That analogy opens up the possi- 
bility to investigate this theory as a constrained system. Having this possibility 
in mind, we propose a Lagrangian to describe this new theory of turbulence 
and, subsequently, analyze it from the symplectic point of view. From this 
analysis, a hidden gauge symmetry is revealed, providing a clear interpretation 
and meaning of the physics behind the metafluid theory. Further, the geomet- 
rical interpretation to the gauge symmetries is discussed and the spectrum for 
3D turbulence computed. 
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I. INTRODUCTION 



The understanding of hydrodynamic turbulence is an important problem for nature sci- 
ence, from both, theoretical and experimental point of view, and has been investigated 
intensively over the last century, but a deep and fully comprehension of the problem 
remains obscure. 

Over the last years, the investigation of turbulent hydrodynamics has experienced a 
revival since turbulence has became a very fruitful research field for theoreticians, that study 
the analogies between turbulence and field theory, critical phenomena and condensed matter 
physics [fi|-|T2|j, renewing the optimism to solve the turbulence problem. 

The dynamics of turbulent viscous fluid is expressed by the Navier-Stokes (NS) equations 
of motion |TJ, which in vectorial form is 

| = ^ xa _vg4) + ,vx (...) 

where u(x, t) is the velocity field, w(x, t) the vorticity field, p(x, t) is the pressure, p the 
density and v the kinematic viscosity. 

The equation of continuity reduces to the requirement that the velocity field is diver- 
genceless for incompressible fluids, i.e, 

V.u = 0, (1.2) 

which are the flows we are interested in this paper. In this context, the hydrodynamic 
turbulence has attracted an enormous interest due to the universal characteristics stressed by 
an incompressible fluid with high Reynolds numbers in the fully developed turbulent regime. 
The Reynolds number, R = LU/v (where L is the integral length-scale of the largest 
eddies and U is a characteristic large-scale velocity), measures the competition between 
convective and diffusive processes in an incompressible fluid described by the NS equations. 
In view of this, the incompressible fluid flow assumes high Reynolds numbers when the 
velocity increases and, consequently, the solution for eqn.( |l.l|) becomes unstable and the 
fluid switches to a new regime of a very complex motion with the velocity varying almost 
randomly and without any noticeable order. To discover the laws describing what exactly 
is going on with the fluid in this turbulent regime is very important to both theoretical and 
applied science. 



Recently, Marmanis [13] has proposed an alternative approach to treat fluid turbulence. 
Based on the analogy between Maxwell electromagnetism and turbulent hydrodynamics, he 
describes the dynamical behavior of average flow quantities in incompressible fluid flows with 
high Reynolds numbers. In this theory, metafluid dynamics, the vorticity (w = V x u) and 
Lamb vector (1 = w x u) are recognized as the kernel of this dynamical theory of turbulence. 



In our paper, we make a further investigation of the physical contents present in that theory, 
in order to furnish a better understanding of turbulence. 

To extract the physical meaning of the metafTuid dynamics, we extend the Marmanis 
analogy in order to propose an appropriate Lagrangian governing the dynamics of incom- 
pressible fluid flows. From this Lagrangian and using the symplectic method ||14|| , the physics 
behind the metafTuid theory is discovered. It is a new and strike result not yet discussed in 
the literature. 

In order to make this work self-consistent, we have organized this paper as follows. In sec- 
tion 2, a brief review of the metafluid dynamics is presented as well as the pertinent physical 
quantities and notation. In section 3, the Lagrangian approach for the metafluid dynamics 
is proposed and some considerations about the dissipation rate of energy is discussed. In 
section 4, the new description for the metafluid dynamics is quantized using the symplectic 
method, which leads to unveil the gauge symmetry. Further, we scrutinize the symmetry of 
the model and show that the gauge invariance is only preserved in some limits, called inertial 
range, and also we identify the physical quantities that are gauge invariant into this region. 
In section 5, the geometric interpretation of the gauge symmetry is given and discussed. In 
section 6, the spectrum of the metafluid theory is computed and the result analyzed and 
compared with the usual results in the literature fl~5 |. In Sec. 7, we stress our conclusion and 
final discussions. We added an appendix with a brief review of the Simplectic formalism. 



II. METAFLUID DYNAMICS 

The problem of turbulence is to find the averaged properties from the solutions of NS 
equations under the constraint of the incompressibility condition, eqn. (|1.2|) , forming a system 
of coupled nonlinear partial differential equations. Due to this, it is a difficult task to get 
a common feature of averaged nonlinear equations, because nonlinearity introduces higher 
order momenta of fluctuation, what is known as closure problem. To overcome this problem, 
Marmanis [|l3l proposed an approximate theory of turbulence based on the analogy between 
Maxwell electromagnetism and turbulent hydrodynamics, where the equations governing 
the dynamic variables become linear and the nonlinearities appear as sources of turbulent 
motion. In this picture, one constructs a system of equations written in terms of the average 
values of the vorticity (w) and Lamb vector (1), instead of the average values of u and p. 
Essentially, it turns the turbulence closure problem into a study of turbulent sources for 
different geometries of interest. 

In this scenario, the study of average quantities of an incompressible fluid at the fully 
developed regime is proposed, given place to a system where the average fields show up in a 
continuum inter-relation and respond as waves to the turbulent sources. To do so, the Lamb 



vector and the vorticity should be taken as the kernel of the turbulent dynamics rather 
through velocity and vorticity fields or velocity and pressure fields. Then whatever parts 
that can not be explicitly expressed as a function of w or 1 only, are gathered and treated as 
source terms. This is done by introducing the concepts of turbulent charge (n) and turbulent 
current (J). The turbulent charge is connected with the Bernoulli energy function, 

through the relation 

n(f,t) = -V 2 $. (2.2) 

In this formalism, the equations of motion describing the behavior of the hydrodynamic 
turbulence are 

V.w = 0, 

— = -V x I + z/V 2 w, 
at 

V.f=n(x,t), (2.3) 
91 

— = c 2 V xw - J(x, t) + uVn(x, t) - z/V 2 /, 

where c 2 = (u 2 ) is the spatial averaged squared velocity, w and I are defined as the averages 
of 1 and w, while the sources J and n are averages of j and n, respectively. Further, the 
Lamb vector can be written in terms of the velocity field and the Bernoulli energy function 
from the NS equations as 

/ = -— -V0 + z/V 2 m, (2.4) 
at 

where u = (u) and = (<&). 

Now, we would like to compare the set of eqn.( [2.3|) with the Maxwell equations with 
sources in vacuum 



16] . Note that these set of equations are identical if B (magnetic field) 
corresponds to w and E (electric field) corresponds to I. Furthermore, the analogy can be 
extended in order to include the potentials as well. In particular, the comparison suggests 
that the vector potential corresponds to the velocity field u and the scalar potential to the 
Bernoulli energy function 0. 

Despite of all resemblance between hydrodynamic turbulence and electromagnetic theory, 
there is a conceptual difference in the identification of the physical entities. In the classical 
electromagnetism, the physical fields are the electric and magnetic fields while the potentials 
are just mathematical artifices. Oppositely, in the metafluid dynamics, the potentials are 
the entities which have physical significance. 



III. LAGRANGIAN APPROACH 



From the geometric point of view and using the Lie algebra, Arnold [T7| showed that 



Euler flow can be described in the Hamiltonian formalism in any dimension. This has a 
lot of interesting consequences for fluid mechanics and has been studied intensively [p~8| , p~9 



However, it is not quite obvious that this process can be used when viscosity is taken in 
account. It is in this scenario where the metafluid dynamics births, revealing a way to find 
a Hamiltonian formalism for a turbulent flow with viscosity. 

In the classical electromagnetism, the Lagrangian density can be written as the difference 
between the square of the electric and magnetic fields, as follows, 

C = \(E 2 -B 2 ). (3.1) 

Using the analogy established between electromagnetism and turbulence, one can write down 
the Lagrangian density for the theory of turbulence as 

£ = ^(f 2 -cW), (3.2) 

that can also be written in terms of velocity field and the Bernoulli energy function, named 
"potentials of theory" , as 




(3.3) 

It is easy to see that this Lagrangian density gives us the equations of motion ( |2.3| ) for the 
homogeneous case (no sources). 

Now, let us consider the case when sources do not vanish. In this case, the sources appear 
in the equations of motion and, as a consequence, an interaction Lagrangian, defined as 

Cint = J-u — ri(p — vu.Vn, (3.4) 

which contains a viscous correction term, is added to the Lagrangian (|3.3| ), providing the 
total Lagrangian density, 

if du \ 2 1 

C = - -V0 - — + z/V 2 u --c 2 (Vxm) 2 + J .u - n<p - vu.Vn. (3.5) 

There is no doubt that the metafluid dynamics can produce some interesting results. Hence, 
let us show that the equations of motion for viscous fluid, given in eqn.( |2.3|) , can be found 
from the Euler-Lagrange equations, where the velocity field (u) and the Bernoulli energy 
function (0) are the canonical variables. From the Lagrangian density, eqn.( |3.5|) , NS equa- 
tions will be obtained in order to demonstrate that the physical contents of turbulence 



in the NS equations description are also present in our construction. The conjugated mo- 
menta of velocity field is computed as 

t) = \ = -El + V0 - v^ 2 u = -l(x, t), (3.6) 

ou(x, t) Ot 



and using (p^), we get the NS equations, namely, 

-u(x, t) = -l(x, t) - V0(f , t) + z/V 2 m(x, t), (3.7) 



dt 

where l(x, t) — w x u. 

One of the simplest consequences of the NS equations ( |1.1|) , modified by the addition of 
the external force /, may be obtained by taking its scalar product with u and integrating 
the result over the space coordinate. The result is the relation 

which expresses the energy balance: the time derivative of energy on the right hand side is 
equal to the difference of the injection rate / f.u and the dissipation of energy rate | / (VtT) 2 . 
Taking averages in the stationary state, we obtain 

(u./) = (^(Vu) 2 )^, (3.9) 

that is the equality of the (intensive) mean injection and the mean dissipation rates of energy. 
The energy injection takes place at the distances of order of the integral scale by induction 
of big scale L eddies. According to the picture of the turbulent flow proposed in 1922 by 



Richardson [20], the big eddies induce still smaller eddies and so on transferring energy 
from large to small distance scales. This process should not lead to a loss of energy until 
sufficiently small distance scales, say, smaller than r\ (dissipative scales of Kolmogorov 0), 
are reached. On scales smaller than rj, the dissipative term z/V 2 u of the NS equation becomes 
important. 

Considering now the intermediate domain of scales called the inertial range: 

n Inertial Range <ti L, 

where the injection rate and the dissipation rate of energy should be negligeable, the La- 
grangian density for the metafluid dynamics becomes 

1 / dil\ 2 1 

C = - f -V0 - — J - -c 2 (V x uf + J.u - ncji. (3.10) 

From this Lagrangian, the analogy between hydrodynamic turbulence and the electro- 
magnetism will be extended. Note that in eqn. (|2.4j ), the fields / and w, as happens with the 



electric and magnetic fields, are components of an antisymmetric second rank tensor, the 
strength tensor, defined as 



which can be written out explicitly for I and w as 



/0 
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The four-vector potential and four- vector current are 

v = (</>, cti), 

r = {en, I), 

respectively. Therefore, one can write (|3.10 ) in term of the F^, J M and V 1 as 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



4 ' c 

It is easily verified through the Euler-Lagrange equations that Q3.14 ) gives the equations 
of motion for the inhomogeneous metafluid dynamics (p. 3D with v = 0. The homogeneous 
equations of motion are obtained from the dual of F^, given by 

1 



a/3- 



(3.15) 



IV. SYMPLECTIC ANALYSIS 

In the last sections, the analogy between turbulent hydrodynamics and Maxwell elec- 
tromagnetism was explored and a gauge field theory to describe the turbulent fluid, called 
metafluid theory, was proposed, suggesting that this theory can be analyzed as a constrained 
system. Both systems are characterized in phase space by the presence of some functions 
that depend on the coordinates and canonical momenta, denominated constraints, which 
restrain the dynamics of the model. There are some methods to handle constrained systems, 
however we analyze the metafluid theory from the symplectic point of view [13]. In the 
appendix we present the symplectic formalism for completeness. 



To implement the symplectic formalism, the Lagrangian density (|3.5| ) will be rewritten 
in its first-order form, given by, 



£(°) = 7 r.£-f/(°), 



(4.1) 



where the canonical momenta and the zeroth-iterative symplectic potential (U^) are 



7T 



+ u — vV u, 



11 - 
U {0) = -re 2 - jr. V0 c 2 (V x m) 2 + z/tt.V 2 m - u.J + 0n + i/u.Vra, (4.2) 

respectively. 

From the set of symplectic variables ^ = (ui, 7Tj, 0) and their respective one- form 
canonical momenta, 



og> = 0, 

4 0) = o, 



(4.3) 



the symplectic matrix is computed as 

/ -Sn 0\ 



/( o) 







5{x-y), 



(4.4) 



V o o o/ 



which is singular, so, has a zero-mode (v^ = (v^, 0, 0)). Contracting this zero-mode with 
the gradient of symplectic potential U^°\ the following constraint appears, 



Q 1 = V.7?(f) +n(x). 



(4.5) 



In agreement with the symplectic formalism, this constraint is introduced into the Lagrangian 
through a Lagrange multiplier A, namely, 



£W = j?M+\n 1 -u<- r > 



(4.6) 



where the symplectic potential density is 

[/« = [/(°) | ni=0 = Itt 2 - ic 2 (V x w) 2 + z/7f.V 2 M - u.J + vu.Vn. (4.7) 

Considering now that the new set of symplectic variables is given in the following order, 
= (ui, 7Tj, A), we have the one-form canonical momenta as, 



(4.8) 



V.7T + n. 



The symplectic matrix is 



( -hj \ 



./ 



(1) 



5ij 







S(x-y), 



(4.9) 



V -d\ J 
that is a singular matrix with a zero-mode given by. 



K : , 0, A 



satisfying the following relation, 



0. 



(4.10) 



(4.11) 



Contracting the zero-mode (v^ 1 *) with the gradient of the symplectic potential (U^), we 
get a new constraint, 

1 



d 3 xvf 



6 



5ui(x, t) 



d 3 y{ — c 2 (V x u) 2 + z/vf.V 2 M + vuSJn - u.J} 
2 



d 3 xvf(x)Ji{x). 



Using ( |4.11| ), fi 2 becomes 



fi 2 = J d 3 xv x (x)diJi(x) = - J d 3 xv x (x)—V 



2 ^x). 



(4.12) 



(4.13) 



At this point, it is important to regard that at the inertial range, despite of the existence of 
viscosity, there is not energy dissipation, then, the Bernoulli energy function ((f)) is constant. 
Due to this, ^2 vanishes and, consequently, the hidden gauge symmetry of the metaffuid 
theory is revealed. 

To finish the symmetry analysis, it is necessary to obtain the infinitesimal gauge trans- 
formation. In the symplectic context, the gauge transformations are generated by the zero- 
mode that does not produce a new constraint. In the present case, the zero-mode given 
by eqn.( |4.1d| ), does not generate a new constraint, consequently, it is the generator of the 
following infinitesimal gauge transformations, 



Sit* = d^, 
8-Ki = 0, 



-e, 



(4.14) 



where — > A and e is an infinitesimal time-dependent parameter. It is easy to verify that 
the Lagrangian (|3.14p is invariant under these transformations, 



5C 



-(<9j7Tj + n)e + (— h + h)e + (d^i + n)i/d 2 e 



0. 



(4.15) 



However, for the Hamiltonian below, 



n 



1 



-7T — 7T. 



we get 



c 2 (V x uf + z/7f.V 2 -u - -u. J + 0n + vu.Vn, 

2 



9 



5ft = e9jJj = -eV 2 -. 

at 



(4.16) 



(4.17) 



We can observe that the gauge invariance is only preserved at the inertial range, where the 
Bernoulli energy function is constant. 

In order to obtain the Dirac brackets among the phase space fields, we have to fix the 
gauge symmetry. It is done introducing a gauge condiction into the kinetical sector of first- 
order Lagrangian through a Lagrange multiplier. So, we choose a gauge fixing term that 
satisfies the condition of incompressibility of fluid, namely, 



X = V.u. 

In view of this, the twice iterated Lagrangian is obtained as 
where the symplectic potential density is 

U(2) = Jjd) 



u 2 



x=0 = -7T + C U.V U + Z/7T.V U — U.J. 



(4.18) 



(4.19) 



(4.20) 



Considering now that the new set of symplectic variables is given in the following order 
d 2) = ( u i, ^i, A, rj), we have 





= 7Ti, 




«S? 


= o, 


(4.21) 


4 2 ' 


= V.7r + n, 






= V.u. 





Then, the corresponding symplectic matrix is obtained as 



/(2 ) = 



/ 










-a? 


\ 




Sij 





-9f 






















V 













/ 



5(x-y). 



(4.22) 



This matrix is nonsingular and, consequently, the corresponding inverse matrix can be deter- 
mined after a straightforward calculation. From the inverse of , the nonvanishing Dirac 
brackets among the phase space fields is automatically identified, namely, 



Sij-^f J5(x-y). (4.23) 

The next step would be the quantization of this constrained theory. Using the well known 
canonical quantization rule ({ , }* — > —i[ , ]), we have, 



[U i {X) ) Uj{y)\ = [TTi(X), 7Tj 

(Qx \ 
Sij-^fUiS-y). (4.24) 

Once we have the canonical quantization rule, we can apply standard Quantum Field 
Theory techniques to find the generating functional and, consequently, the correlation func- 



tions and all physical quantities []2T| , [22| one wish. 



V. GEOMETRIC ASPECT OF THE GAUGE SYMMETRY 

In this section, we give a geometrical interpretation to the gauge symmetry present on 
the metafluid theory (gauge field theory). The description of fundamental particle interac- 
tion, with the assistance of the gauge field theory, introduces extra degrees of freedom into 
the theory, that manifests itself in the singular nature of the respective Lagrangian or the 
presence of first class constraint in the equivalent Hamiltonian formulation. In this case, 
the phase space is larger than the physical one, which is a hypersurface determined by the 
constraints of the theory. In this gauge invariant scenario, the gauge potentials form an 
overcomplete basis and the gauge fields, which can be connected by an infinitesimal trans- 
formation, describe the same physical state. Thus, the gauge potentials are separated into 
equivalence classes with respect to the gauge group action, where each one denotes an orbit 
in the gauge field configuration. Transitions along(vertical) the orbits correspond to pure 
gauge transformations, consequently, these paths have no physical significance. Oppositely, 
perpendicular(horizontal) paths to the orbits describe the time evolution of the physical sys- 
tem, then, they are physically important to the theory. Hence, to find out the equations of 
motion for the physical fields consist in solving the problem of constructing the horizontal 
paths, which is computed just doing a correct definition of the metric in the orbit space. To 



do so, we follow process developed in ||23|| , where the authors demonstrated that the physical 
(orbit) space is equipped with a natural projective metric. 

In view of this, our task is to get a singular Lagrangian, dynamically equivalent to ( |3.14j ), 



to govern the metafluid theory. It is achieved eliminating the Bernoulli energy function (0) 
from the Lagrangian (|3. 14 ) . It can be done since this field has no dynamics, as demonstrated 



in the last section. To eliminate this field from ( |3.14j ), we use the Euler-Lagrange equation 
of motion for <ft, given by, 



i/->\ diiii + n 

4>{x) = ^2 h vdiUi. (5.1) 



Bringing back this result into the Lagrangian ( 3.14 ) and rewriting it in terms of field com- 
ponents, we get 



C = -UiMijiij + (diiii)— + 2Q2~ v(d 2 Ui)M, 



LljUj 



v 2 1 
+ —(d^Miju, + -c 2 (d 2 Ui )M ijUj + mJi, (5.2) 

where the metric of configuration space M^, reads as 

M l3 = * y - ^, (5.3) 

being a singular matrix which has d{ as eigenvectors with zero eigenvalues, 

diMij = 0. (5.4) 

As the gauge orbits lay down by the eigenvectors, are vertical to the physical hypersurface 
defined by the projective metric Mij(M 2 = M), the infinitesimal gauge transformation for 
the velocity potential can be computed as 

5ui = die, (5.5) 

which was also obtained in symplectic context. 

At this stage, we would like to make some comments about the geometry associated to 
the gauge symmetry. From the Lagrangian ( |5.2|) , we identify a singular projective metric 
which defines a physical surface (orbit space) of the metafluid theory. Since this metric does 
not depend on the phase space fields, the curvature tensor is null, therefore, the orbit space 
is flat. 

Now, we investigate the metafluid theory described by the Lagrangian ( |5.2p using the 
projector method [pi ]. From eqn.( |5.2| ), the canonically conjugated momenta to the field Ui 
are computed as 

iii = Mijiij _ "qV _ vMijd 2 Uy (5.6) 

Contracting these canonical momenta with the eigenvectors <9,, the primary constraint is 
determined as 

X = di-Ki + n (5.7) 



which agrees with the results obtained in the symplectic analysis. 

From the Lagrangian (|5.2| ), the Euler-Lagrange vector is obtained as 



E k = M ki Uj - c 2 M kj d 2 Uj -Jk-^r- v 2 M kj 8\. (5.8) 

Using the projector method [p4j| , the equations of motion that satisfy the constraints are 
projected by the singular metric, namely, 

M lk E k = 0, (5.9) 

and, then, we get the wave equation for the velocity [O, given by, 



Mijiij - c 2 Mij^Uj - MijJj - v 2 M i ^u j = 0, 

ui -c 2 d 2 ui - -u 2 8 4 ui = ^ (5.10) 

where the velocity and current were redefined as uf = M^Uj and J±- = M^Jj, respectively, 
and _L denotes the transverse fields (gauge invariant fields). This set of equations together 
with the constraint condition, eqn. (|5.7| ), define the metafluid dynamics with transverse fields. 
Due to this, the incompressibility condition is satisfied automatically, namely, 

diUi~ = diMijUj 

d i m ± = 0- (5.H) 

Contracting the eigenvector <9j with the Euler-Lagrange vector, 

d k E k = d k {M kj u 3 ) - c 2 d k (M k ,d\) - d k J k - d k ^ - u 2 d k (M kj d\), (5.12) 
we get the continuity equation, given by, 

d k E k =n + d k J k = 0. (5.13) 

As n = — V 2 0, we obtain 

- V 2 + d k J k = 0, 

d k J k = V 2 = 0, (5.14) 

because there is no energy dissipation in the metafluid theory, since it was constructed at the 
inertial range. Therefore, at this region, the current satisfies the divergenceless condition. 



VI. SPECTRUM OF THE THEORY 

The close resemblance of the equations (|2.3|) to the equations of electromagnetism im- 
mediately suggests that turbulent flows have a hidden wave character in them. Indeed, 
substituting ( |2.4j ) and w = V x u in the last equations of (|2.3|) , it can be shown that the 
average velocity field obeys a wave equation, given by 



— = c 2 V 2 u + z/ 2 VV+ J tr , (6.1) 

— # 

where J tr is the turbulent transverse current density. 

In the homogeneous turbulence, the turbulent sources do vanish and eqn.( |6.1| ) has the 
following form 

?f = cV« + . 2 V 4 « (6.2) 
ot 

or 

Vu = 6 (6.3) 

where the T> operator is given by 

V = d 2 t - c 2 V 2 - u 2 V\ (6.4) 
Now, we can get the Feynmann propagator or the Green function for eqn. (|6.3|) , 

iSf^x-x') = (0\T( Ul (x) Uj (x'))\0}, 

^[x x)- J ^ kl _^ ^ % J, (6.5) 



with 



T{ui(x)uj(x')) = 8(t — t')ui(x)uj(x) + 9(t' — t)uj(x')ui(x), 
k = ko — k, 

x = xo — x, (6.6) 

2 2 f 2 2 f 4 

uj h = c k —vk, 



where the last equation is the dispersion relation to the wave equation (|6.1| ) Hl3 |. 

The calculation of S^(x — a;') can be performed as usual and the result for equal time 
Feynmann propagator is given by 

Now, summing over for i = j, we get 

3 _ /■ H 3 ? p ik.(x-x') 

= / (^3—- (6-8) 

Therefore, eqn.( |6.8|) can be identified with the two-point correlation function of the ve- 
locity field at equal time. So 



(u(x).u(0)) = ^iSf^x-x') = / —- (6.9) 

or in the Fourier space, we get 

(u{k).u{-k)) = J d 3 x e~ ils (u(x).u(0)) = u^ 1 . (6.10) 

With this calculation, we showed (u(k).u(—k)) oc E{k) oc k~ x . Such result would be 
expected because of the analogy between electromagnetism and the metafluid dynamics. The 
reason why in this model we get this exponent instead of the value —5/3 found in literature 
is that the linearization implies the decomposition of the field into various oscillators do not 
invoke any interaction among them. To get the exponent found in the literature on can 



proceede as done by Zakharov et al [^,[L3|], where a perturbation theory was performed. 
Despite of the value obtained in this formalism, one can see the dispertion relation u(k) 
presents the physics expected for a turbulent flow once, u(k) can be real, with positive or 
negative, or even imaginary. This means in the context of Fourier space that one can have 
waves with positive or negative energy and that u(k) been imaginary gives a dissipative 
behavior, but that can happen only for sufficiently high wavenumbers. 



VII. CONCLUSION 

In this paper, we proposed to investigate incompressible turbulent hydrodynamics in the 
context of metafluid dynamics in order to open up the possibility to apply all the machinery 
very well known in quantum field theory. This investigation was possible due to the analogy 
of incompressible turbulent hydrodynamics with Maxwell electromagnetism. Exploiting this 
analogy, the Dirac brackets was computed through the symplectic formalism, which revealed 
a hidden symmetry presents on the metafluid gauge theory, only preserved in inertial range. 
Afterward, the geometrical meaning of the gauge symmetry was given, showing that the orbit 
space is flat and, consequently, the divergenceless condition (V.J 1 - = 0) is naturally obtained, 
at least in the inertial range. Subsequently, the wave equations for the transverse (horizontal) 
velocity field were computed, which was also obtained in |Tj|, but written in terms of the 
velocity field, both horizontal (gauge invariant) and vertical (gauge variant) components. To 
finish, the 3D spectrum of metafluid theory was computed, leading to a proper result which 
does not agree with the usual one computed in the incompressible turbulent hydrodynamics 
context ||15|| . It seems that the approximate approach delineated in |l3j and extended here 



with the introduction of the Lagrangian, given in eqn. fl3.14| ), change the spectrum. We 



suspect that some singularities due to the nonlinear nature of the incompressible turbulent 
hydrodynamics are eliminated by the analogy process used to linearize the theory. We believe 



that the contribution of these singularities in the metafluid theory can appear if the orbit 
space has a curvature. Despite of not getting the usual exponent for the energy spectrum, 
we showed that with this formalism, one gets the right behavior for the turbulent flow, in 
the inertial range. 
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APPENDIX A: SYMPLECTIC FORMALISM 



Faddeev and Jackiw Eg] and Barcelos and Wotzasek [27] showed how to implement the 



constraints directly into the canonical part of the first-order Lagrangian once this method is 
applied to first order Lagrangians. In order to systematize the symplectic gauge formalism, 
a general noninvariant mechanical model that has its dynamics governed by a Lagrangian 
£(a,i,a,i, t) (with i = 1,2, ...,N) is considered, where and dj are the space and velocities 
variables respectively. Notice that this consideration does not lead to lost of generality. In 
the symplectic method, the first-order Lagrangian written in terms of the sympletic variables 
£W(a,i,Pi) (with a = 1, 2, 2N), is required 

£(°) = 4°)ei 0) -[/ (0) , (ai) 

where is the one-form canonical momenta, (0) indicates that it is the zeroth-iterative 
Lagrangian and, U^°\ the sympletic potencial. After that, the sympletic tensor, defined as 

,<o> _ dJl _ &4£) 

is computed. Since this sympletic matrix is singular, it has a zero-mode (z/°)) that generates 
a new constraint when contracted with the gradient of potencial, namely, 

n<°> - ^gf ■ («) 

Through a Lagrange multiplier A, this constraint is introduced into the zeroth-iterative 
Lagrangian ( |A1| ), generating the next one 

£(i) = 4°)^ 0) -t/ (0) + A^ ), 

= Ag>g?> - (A4) 



where 



A y = 4°) + A 



fi(°)=0> 

^M^U), (A5) 



The first-iterative sympletic tensor is computed as 

Since this tensor is nonsingular, the iterative process stops and the Diras brackets among 
the phase-space variables are obtained from the inverse matrix. On the other hand, if the 
tensor is singular, a new constraint arises and the iterative process goes on. 
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